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The CP-restoring phase transition at = tt and high temperature is investigated using two related 
models that aim to describe the low-energy phenomenology of QCD, the NJL model and the linear 
sigma model coupled to quarks. Despite many similarities between the models, different predictions 
for the order of the phase transition result. Using the Landau- Ginzburg formalism, the origin of this 
difference is traced back to a non-analytic vacuum term at zero temperature that is present in the 
NJL model, but usually not included in the linear sigma model. Due to the absence of explicit CP 
violation, this term always alters the qualitative aspects of the high temperature phase transition 
at 9 = n, just as for 6* = in the chiral limit. 

PACS numbers: 12.39.-x,11.30.Er,11.30.Rd 

It is well known that there is a possibility of CP violation in the strong interaction due to instanton contributions. 

n 2 

These contributions are incorporated in the QCD Lagrangian through the topological F F-teiui, where 6 is the 
QCD vacuum angle. This term violates CP, unless = mod tt. The case 6* = tt is special, because then Dashen's 
phenomenon can occur, i.e., spontaneous CP violation at = tt [ijl. 
I—, ^ From experiments it is known that in nature 9 is very small [2, 0, 0, [H| ■ The reason for this is unknown and is 
commonly referred to as the strong CP problem. However, it has been argued that in heavy-ion collisions meta-stable 
Qh^ CP-violating states could be created corresponding to states with an effective 9 [1, 0, H, H, [l^j [HI [3 ■ Studying 
O • the behavior of the strong interactions at nonzero 9 is therefore of interest and has been done quite extensively using 
<]J ' chiral Lagrangians, see for example Refs pi. [Tsl. [TgI. [TtI. [Tsl. [l9l [20t . 

, Recently the 6'-dependence of two models describing the chiral dynamics of low energy QCD have been studied, the 
NJL model [2lj and the linear sigma model coupled to quarks (LSMg) . In both models the effects of instantons 
^f-^ are included through an additional interaction, the 't Hooft determinant interaction (23l. [23|. It was found that both 
^ ' models exhibit Dashen's phenomenon, which turns out to be temperature dependent. This is to be expected, since at 
\ high temperature the effects of instantons, which are needed for the CP violation, are exponentially suppressed [2^ . 
^^O . In both models the spontaneous CP violation at 6* = tt disappears at a critical temperature between 100 and 200 
^2 , MeV, however, the order of the phase transition differs. In case of the NJL model the transition is of second order, 
. whereas in the LSMg' model it is of first order. Clearly this difference is important, because a first order transition 
■ allows meta-stable phases, in contrast to a second order transition. 

, Although the NJL and LSMg model are not the same, they are closely related. Eguchi [2^ has shown that when the 
NJL model is bosonized, a linear sigma model is obtained (see also [23|). However, the effects of quarks are treated 
differently in both models, which was already discussed in Ref. [1^ for 9 = 0. In the case of the LSMg model the effects 
of the quarks are usually only taken into account for nonzero temperatures, whereas in the NJL model their effects 
are necessarily incorporated also at zero temperature. Ref. ^285 found that the order of the chiral symmetry restoring 
phase transition ad, 9 = was the same in both models, but the critical temperatures differ. While the qualitative 
' aspects of the phase transition are similar at 6* = 0, this is not the case for the high temperature CP-restoring phase 
transition at 9 = n as we will discuss in detail. We should mention here that the situation at 9 = depends on the 
amount of explicit chiral symmetry breaking. In Ref. [29] it was observed that when the pion mass is reduced in order 
to study the chiral limit, neglecting the effects of the quarks at zero temperature can affect the order of the high 
temperature phase transition at 9 = too. 

Although there is a CP-restoring phase transition at high chemical potential also, in this paper we will restrict to 
the temperature dependence of this phase transition at 9 = n, because there the differences between the two models 
are most pronounced. The paper is organized as follows. First, the effective potentials of both models are analyzed 
analytically, which will allow the determination of the order of the phase transitions using standard Landau-Ginzburg 
type of arguments. A comparison to numerical results obtained earlier corroborates these conclusions. Subsequently, 
we will discuss the bosonification procedure of Eguchi, which relates the NJL model to a linear sigma model and 
allows us to further pinpoint the origin of the similarities and differences with the LSMg model. 



o 



X 



2 



I. NJL MODEL 

The Nambu-Jona-Lasinio (NJL) model, introduced in Refs. [13, HH, is a model for low energy QCD that contains 
four-point interactions between the quarks. In this paper the following form of the NJL model is used, in the notation 
of Ref . [U 

Cnjl = -il^ (il'^d^ ~ m)Tp + Cqq + Cdct, (1) 

where m is the current quark mass. In contrast to Ref. 21], here the up and down quark masses are taken equal, 
which matters little for our present purposes. Furthermore, 

Cqq = Gl [(^TaV-)' + (V5raZ75^)'] , (2) 

is the attractive part of the qq channel of the Fierz transformed color current-current interaction ^32l| and 

-Cdct = 8G2e''det(Vifl^L) +h.c., (3) 

is the 't Hooft determinant interaction which depends on the QCD vacuum angle 9 and describes the effects of 
instantons [H, [131 ■ In the literature Gi and G2 are often taken equal, which at 6* = means that the low energy 
spectrum consists of a and tt fields only, but here we will allow them to be different. We will restrict to the two 
flavor case, using with a = 0, ...,3 as generators of U(2). We will not consider nonzero baryon or isospin chemical 
potential. 

The symmetry structure of the NJL model is very similar to that of QCD. In the absence of quark masses 
and the instanton interaction, there is a global SU(3)cXU(2)ixU(2)i^-symmetry. The instanton interaction 
breaks it to SU(3)cXSU(2)ixSU(2)ij;xU(l)_B. For nonzero, but equal quark masses this symmetry is reduced to 

SU(3)cXSU(2)yXU(l)i3. 

We choose the parameters the same way as in Refs. [2l|, [s^]. This means we write 

Gl = (1 - c)Go, G2 = cGo, (4) 

where the parameter c controls the instanton interaction, while the value for the quark condensate at = (which 
is determined by the combination Gi -\- G2) is kept fixed. For our numerical studies we will use the following values 
for the parameters: m — 6 MeV, a three-dimensional momentum UV cut-off A = 590 MeV/c and GqA^ = 2.435. 
These values lead to a pion mass of 140.2 MeV, a pion decay constant of 92.6 MeV and finally, a quark condensate 
(uu) — (dd) — (—241.5 MeV)'^ [3^, all in reasonable agreement with experimental determinations. 



A. The effective potential 



To calculate the ground state of the theory, the effective potential has to be minimized. In this section the effective 
potential is calculated in the mean-field approximation. In the following we will only consider the case of unbroken 
isospin symmetry, such that only nonzero (V'V') and/or (ipi^/c,^) can arise. At = only (ip'ip) becomes nonzero. A 
nonzero (i/jij^tp) signals that CP invariance is broken, i.e., it serves as an order parameter for the CP-violating phase. 

To obtain the effective potential in the mean-field approximation, first the interaction terms are "linearized" in 
the presence of the ('i/'^A) and {ipi^^ili) condensates (this is equivalent to the procedure with a Hubbard-Stratonovich 
transformation used in Ref. [21| ) 

(^V)' - 2 (V^V-) - (^V-)" , 
{ip'ip){'ipij5ip) ~ {ipip) ipijii' + ii'nsip) i'lp - {H-') (V'nsV') , (5) 

leading to 

/•vac ^J,(^^,f^» i,A\j, (Gl -G2 cos6>)q^ (Gl -KGz cos6>)/3g (G2 sin6>)Qo/3o 
^NJL d^-M)^ 4(G? - Gl) 4(G2 - G^) 2(G2 - G^) ' 



where Ai = {ra + a^) + /3o*75 and 



ao = -2(Gi-t-G2Cos6i)(V5V'> + 2G2sin6'(V;i75V'), 

/3o = -2(Gi-G2Cos0)(V;i75V'> + 2G2sin0(V;^). (7) 
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This Lagrangian is quadratic in the quark fields, so the integration can be performed. After going to imaginary time 
the thermal effective potential in the mean-field approximation is obtained [33| 



q;§(Gi - G2 cos6') /3^(Gi + G2COS6I) G2ao/9osin( 
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po = (2n+l)7rT 



d^P 
(27r)3 



log det K, 



(9) 



and where K is the inverse quark propagator, 



K = {ijopo + -fiPi) - M. 



(10) 



In order to calculate the effective potential, it is convenient to multiply K with 70, which does not change the 
determinant, but gives a new matrix K with ip^s on the diagonal. It follows that deiK — Y^^^l {Xi — ipo), where Xi 
are the eigenvalues of K with po = 0. Because of the symmetries of the inverse propagator, half of the eigenvalues are 



equal to Ep = y^p^ + i\JP and the other half to Ep 
over the Matsubara frequencies, we obtain 



+ M"^, with = (m + ao)^ +/3o- After the summation 



(11) 



At T = this integral can be performed analytically. A conventional non-covariant three-dimensional UV cut-off is 
used to regularize the integral and yields: 



\M\ (MHogff^ 
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where the degeneracy factor i/q = 2A. 



B. The CP-restoring phase transition 

In this section the high-T CP-restoring phase transition at = tt is investigated. As was shown in Ref. [2l| the 
phenomenon of spontaneous CP violation is governed by the stren gth c of the 't Hooft determinant interaction. It will 
be assumed that c is 0.2, which following the arguments of Ref. [33| is considered realistic. But in fact, the critical 
temperature is too very good approximation c-independent for c above 0.05, as can be seen from the (T, c) phase 
diagram given in Ref. [21] . 

We will start with a numerical minimization as a function of the temperature, the results of which, together with 
those for the LSMg model, are shown in Fig. [TJ One observes that the critical temperature of the NJL model is 
significantly larger than the one of the linear sigma model, in agreement with the results of Ref. [2^ for the chiral 
phase transition at = 0. Furthermore, the order of the phase transition is clearly different, contrary to the results 
of Ref. m for = 0. 

Next we will derive an analytic expression for the effective potential for the NJL model. Two important observations 
which can be made from the numerical study will be help. First, we note that ao is very small and constant as long 
as Po is nonzero, which allows us to approximate ^ /3g. Furthermore, /3o and hence M can be considered much 
smaller than ttT and A, allowing expansions. These observations simplify our study considerably. 

The phase transition occurs for M much smaller than A, so Eq. can be expanded in M/A at T = 0: 

M-^ logM^ MMog(4A2) A^AP A"^ 1 

^ 1287r2 " Wn^ ^ 16^ ^ " ' ' ^ ' 

For the phase transition, the non-analytic term logM^ turns out to be very relevant. We will see that it is exactly 
the absence of this term at finite temperatures in the NJL model that causes the differences between the two models. 

Usually the temperature-dependent part of the potential has to be evaluated numerically, however when M < ttT 
the integral can be expanded in M/T. As can be inferred from Fig. [T] it is exactly this regime which is relevant for the 
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FIG. 1: The temperature dependence of the condensates in the NJL and linear sigma model. 



phase transition. Note that the temperature-dependent part of the potential is UV finite, which means that for this 
part the cut-off can be taken to infinity. In Ref. [2l| this was not done, leading to a slightly larger critical temperature. 
Performing the expansion, we obtain |35| 



d^p 
(2^ 



Tlog 1 



720 



48 



32^ 



3 1 M2 , 

7B - ^ + 2 log ^ " log 



(14) 



From this expansion one can see that also the temperature dependence contains a logarithmic term, that will precisely 
cancel the one of Eq. when added together. 



Using that « /3^, we end up with the effective potential 



where 



(Ttt^T^ + 45A4) 



^njl(T) = - 
Snjl(T) - 

CNJL(r) = 



7207r2 
48^1^ 

(log (4A^) -logT^) 
647r2 



1 

4g;^' 

(-1 + 7g - log7r)i^g 
327r2 



(15) 

(16) 
(17) 
(18) 



One observes that the logarithm at zero temperature is cancelled by the logarithm in the temperature dependence. 
As long as /3o < ttT, A the potential contains no logarithms and is fully analytic. We note that this expression is the 
same as the chiral limit at = 0, with /3o replaced by olq. 

The phase transition occurs when -Bnjl(T') changes sign. As the potential is symmetric and quartic in the order 
parameter, we conclude (following Landau- Ginzburg arguments) that the phase transition is of second order, which 
the numerical analysis corroborates. The critical temperature is equal to 



nNJL 



/3i/gGoA2 - 127r2 



185 MeV. 



(19) 



As long as T < ^ exp(— 1 -I- 7^;) = 246 MeV, Cnjl is positive, such that higher order terms in /3o are not needed in 
the analysis. 
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II. LSMg MODEL 

The linear sigma model coupled to quarks, like the NJL model, is an effective low-energy model for QCD (28l[29l. [36l. 
[37I , similar in form to the Gell-Mann-Levy model [38;] . It is a hybrid model that includes both meson and constituent 
quark degrees of freedom, the latter only at nonzero temperature however. As was the case in the NJL model, the 
effects of instantons are included via the 't Hooft determinant interaction. In this paper the analysis of Ref. [2^ is 
followed. 

We will start with the T = case, when only mesons are considered. The Lagrangian, which contains all Lorentz 
invariant terms allowed by symmetry and renormalizability has the following form, using a slightly different notation 
than Ref. 

Cls = iTr(a^0taM^) + ^Tr(0t<^) _ ^[Tr(0t0)]2 _ ^Tr[(0t<^)2] + |[e^Vet(0) + e-^Vet(0t)] 

+ iTr[^(<^ + </.t)], (20) 
where (f> is chiral field, defined as 

cl) = + If]) + ^{ao + itt) ■ T . (21) 

The Lagrangian incorporates both spontaneous and explicit breaking of chiral symmetry, the latter through the term 
proportional to H. To study this symmetry breaking, we can concentrate on the potential corresponding to Eq. (j20p . 
expressed in the meson fields 

-■^ cos6 (o-^ + TT^ - 77^ - al) 
+K sin9 [arj — tt • ao) — Ha 

2A 

+ ^(o-ao + ?77r + 7r X ao)^ . (22) 

The spontaneous symmetry breaking manifests itself through nonzero a and 77 condensates and are obtained by 
minimizing the potential. We allow for these condensates by shifting the fields 

cr ^ (To + .5, V ^ Vo + h, (23) 

where ao and rjo are the values that minimize the potential and s and h are the fluctuations. These ctq and 770 are 
proportional to the condensates ao and Pq of the NJL model, respectively. 

The potential can now be split in two parts, a vacuum part and one that depends on the fluctuations, i.e., 

vl^"-K7'^=' + v^r- (24) 

First we concentrate on the vacuum part, which is given by the following expression: 

T,vac,T=0 2 2\2 rj , ^/ 2 2\2 , ■ n ,'^22 '^/4|4\ /r,r\ 

^LS =^i<^o-'^e) - Hao + -ivo -ug) + k sm0 aoVo + ^ctoVq - ^(ve + Ug), (25) 
where we have defined the combination of couplings A = Ai + A2/2, and follow the notation of Ref. [23 |: 

2 ^'^ + K CO80 2 2 2k 

vl = ^ ; ul = v^g-—cos9 . (26) 

This part of the potential determines the phase structure and has to be compared with the NJL expression (jl3p . The 
main difference is that this potential is fully analytic and does not contain any logarithmic terms. 
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The part of the potential that depends on the fluctuations is used to determine the parameters /i^, k, iJ, Ai and 
A2 in Eq. (|20p . They are obtained by fitting the masses contained in V^^'^ and the pion decay constant at 6* = such 
that the model reproduces the low-energy phenomenology of QCD. At 6* = Vls*^ has the following form 



Vfluc 
LS 



r22i 22i 22, 2 21 

[m^TV +m^s +m^ri +m„^aoJ 



Ai + -A2 



aos Is 



Ai + ^A2 ) (Josal + X'lCfQrjTT ■ Oq 



iA2 

2 



(sao 



7/7r)^ + (tt X aof 



(27) 



The masses depend on the parameters of the model as follows: 



ao 
™2 



A2)cr, 



A2)cro' 



(Ai 



^A2)a^ 



(28) 



The mass values used are: = 138MeV, = 600MeV, niaa = 980MeV, and = 574MeV. 

At nonzero 9, rjo becomes nonzero, which alters the mass relations. Furthermore, cross terms like ar] become nonzero, 
signalling that the mass eigenstates are no longer CP eigenstates, as discussed for the NJL model in Ref. [2l|. As 
a consequence, the cr-field mixes with the ?7-field and the 7r-field mixes with the ao-field. We will not give these 
expressions explicitly here. 



A. Nonzero temperature 



In the LSMg model the quarks start to contribute at nonzero temperatures. In fact, it is assumed that all the 
temperature dependence comes from the quarks. In Ref. [2^ it is argued that this approach is more justified for 
studying high T phenomena than considering only thermal fluctuations of the meson fields, because at high T con- 
stituent quarks become light and mesonic excitations heavy. For the study of the chiral phase transition at 6 = this 
approach yields results that are qualitatively similar to those of the NJL model. 

The part of the LSMg Lagrangian that depends on the quark fields is: 

Cq = tp [i<^ - g{(T + ij^r-j + ao ■ t + ij^iv ■ r)] ■tp. (29) 

The quark thermal fiuctuations are incorporated in the effective potential for the mesonic sector, by means of inte- 
grating out the quarks to one loop [23| . The resulting quark contribution to the potential is given by 

Vj = -«/(g,riog(l + e-^/-). (30) 

This expression is equal to the temperature dependent part of the NJL model Eq. (fH|) , with again Ep = y^p'^ + M"^ 
and the constituent quark mass M depends on the vacuum expectation values of the meson fields in the following 
way: M = SVv'o^^o)' where g is the Yukawa coupling between the quarks and the mesons. A reasonable value for 
the constituent quark mass aX 9 — Q fixes this coupling constant. In Ref. [23] (and here) g = 3.3 is used, which leads 
to a cross-over for the chiral phase transition as a function of temperature at = and to a constituent quark mass 
of approximately 1/3 of the nuclcon mass. 



B. The phase transition 



With all parameters fixed, we can study the CP-resto ring phase transition at ^? = tt in the LSMg model. This was 
studied in detail, along with other values for 9, in Ref. [22|. There also the effect of a magnetic field was discussed, 
which we will not take into account. 
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We are now going to follow the same procedure as for the NJL model to study the details of the phase transition. 
Again, we start the discussion with numerical results of the minimization of the effective potential, this time the results 
of Ref. [l^l- They are shown in Fig. [TJ From this figure, two simplifying assumptions can be inferred. First, as was 
the case for the NJL model, in the neighborhood of the phase transition M < ttT, allowing Eq. (|30p to be expanded 
in M/T as in Eq. (|14[) . Second, ctq is much smaller than -qo which means that we can neglect the ao-dependence. This 
assumption leads to a small error near 770 ~ 0, but as we checked explicitly this is not important since the structure 
of the extrema of the potential is not altered. 

Summing the contributions at zero and nonzero temperature gives the following form for the effective potential 

VltiT) = A^siT) + BUT)ril + Cls(T)%' + ^ls^^ log%', (31) 

where 
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A^T) = -—n^T\, (32) 
1 

48 



720 ^ 

Bls{T) = j-{9'T',^g-24(p' + K)), (33) 



The form of this potential is clearly different from the one of the NJL model Eq. the difference being the 

uncanceled logarithmic term. This term proportional to Uls will always cause the phase transition to be of first 
order. As observed for the NJL model, also in this case the potential is exactly the same as the chiral limit at 6* = 0, 
with ?7o replaced by (Tq. Beyond the chiral limit the explicit symmetry breaking term ^ Hao (which has no analogue 
at = tt) will change the first order transition into a cross-over, unless the Yukawa coupling g is increased sufficiently 
[29I IstI . We conclude that the absence of explicit CP violation through a linear term in 770 at 6* = tt lies at the heart 
of the difference between the observations made here and those in Ref. 28] . 

Like in the NJL model, it is the sign flip of -Bls that modifies the structure of the minima. But instead of a phase 
transition, now a meta-stable state develops at rjo — 0. When i?Ls(r) becomes larger than 2DLsexp(— | — '"dlJ^ ) 
the original minimum disappears. Between the two spinodals the minimum jumps, signalling a first order transition. 

When the parameters of Ref. [l^ are used, we obtain the following values for the spinodals: 118 MeV and 129 
MeV. To find the exact point of the phase transition, the potential has to be minimized numerically, giving a critical 
temperature of 126.4 MeV. As already noted, this is significantly lower than T^'^^, but the specific values depend on 
the parameter choices made. As should be clear from the previous discussion, choosing different parameters would not 
affect the conclusion about the different orders of the phase transition, at least as long as M < ttT, A and k > — /i^ 
(equivalently, > 3m^ at T = 0). 



III. RELATION BETWEEN THE NJL AND LSM<j MODEL 

As mentioned, the LSMg model is a hybrid model for mesons, which are coupled to quarks at nonzero temperature, 
and the NJL model is a quark model, where the bosonic states of quark-antiquark fields are interpreted as mesons. 
Eguchi [26[ has shown how to derive from the Lagrangian of the NJL model a Lagrangian for the mesonic excitations 
for G2 = 0- This bosonification procedure is reviewed in Ref. (27j . Here the corresponding meson Lagrangian will be 
derived for G2 7^ 0, which was also studied in Ref. [s^ in the chiral limit. 

The situation will be reviewed for — 0, when only the tjjip receives a vacuum expectation. We start with the 
Lagrangian given in Eq. ([T]). The generating functional is given by the standard expression 

^] = ^ / ^^^^ exp (^i J d^x [£njl(^, ^) + + ) , (36) 

where ^ and ^ are the antifermion and fermion sources and A'^ is a normalization factor which will be suppressed from 
now on. Next we introduce auxiliary fields ct, 77, tt and Oq and a new Lagrangian C such that the effective potential 
can be written as 

i]^ j VtpV^pVaVjjVirVao exp (^i J d^x [^njlCV^, V') + V-C + ^V'] j , (37) 
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with 



and 



-^NJL =il}[i^ -m- g{a + i'ynv + ao • r + hstt • t)] - (cr^ + tt^) - ^5/i2 (??^ + ag) , (38) 



2(Gi+G2)' 



2(Gi - G2) 



(39) 



Here 5 is again the Yukawa coupling between the quarks and mesons, which in the case of the NJL model can be 
evaluated. It is equal to 



5-2 = _47V,z 



1 



(27r)4 (p2 _ M2)2 



(40) 



which requires some rcgularization. 

Integrating out the quarks gives the following generating functional 



VaVriD-KVaQ exp I IiSnjl + W d^a;^ 



1 



m- g{8 + ijsT] + ao • T + i^5-!v ■ r) 



(41) 



where the action iSnjl is equal to 
1 



'NJL 



--5(4 {a^ + TT^) - -Si4 {rf + al) 



— iTr log [i^ — m — g{a + i'y^r] + ao ■ t + 17577 • r)] . (42) 



Assuming that only the cr-field receives a vacuum expectation value ctq, i.e., a = ao + s, the action can be split into a 
vacuum part and a part that depends on the fluctuations, which are the mesons s, r], tt, oq: 



cvac I ofluc 

NJL — «->NJL "T '-'NJL) 



(43) 



with 



cvac 
^NJL 



:»fluc 



-I 



•^NJL 



d'3 



-iTrlog [i^-M], 
1 



-Sill {s^ + + TT^) - -SiA {rf + af) 



-iTr log 



1 



id- M 



g{s + ij^t] + ao • T + ij^n ■ r) 



(44) 



and the constituent quark mass M = m + gao. In order to obtain a local action for the meson fields, the nonlocal 
fermionic determinant in «S^jl is rewritten using a derivative expansion: 



— iTi log 



I {s + ij5V + ao • r + 1757? • r) 



(45) 



n=l 



where 



[/(") = Itt 



id-M 



g{s + «75?7 + ao • T + ^7577 • r) 



(46) 



From power counting we note that [/("^ with n > 5 arc convergent and the rest is divergent. Evaluating and retaining 
only the divergent parts of the /7^"^ with n = 1,2,3,4 we end up with the following Lagrangian, which integrated 
over all space yields iS^jl- 

= I [{d.sf + {d,nf + {d.alf + {d,^f] \ [mW + mis' + + m^„a^] 



-93S (s^ + TT^ + r?^ + 3ao) - 2gi'q-K ■ ao - ^54 (s^ + tt^ + r/^ + ao)^ 



-254 



(stto + r]7r)' + (tt X ao)^ 



(47) 
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The masses and coupling constants have the following values 



1 _ /o ml 



lYl ^ — ^ 

2Go/o ^0 M 2GoIo 

2 1 ^2 

"^^ 2(l-2c)Go/o lo' 



2M 



53 = = 2M54 

-^0 



1/2 



94 



1 

io' 



7-1/2 

where /g and /2 are two divergent integrals, equal to 



9 = (48) 



d^p 1 

(2^ (p2 - Af 2)2 ' 



If these integrals are regularized using the three-dimensional UV cut-off, the resulting masses are equal to the ones 
obtained using the random phase approximation used in Ref. (2T| (where the dependence on the external momentum 
of the generalized Iq defined in, for example, Ref. [131 has been neglected). The Lagrangian P7|) without the Uq and 
77- fields was also given in Ref. [4^ . In the chiral limit the results agree with those of Ref. 39] . 

Eq. (|47)) is equal to the fluctuation part of the linear sigma model Lagrangian ((20| using the following parameters 

Ai = 0, 
A2 = 4//0, 
..2 - 2M^ ^ 



2(1 - 2c)Go/o 

c 



2(l-2c)Go/o' 
m 



H - — . (50) 

2G0/0/' 

Although the bosonification of the NJL model yields Ai = 0, this is of no consequence for the order of the phase 
transition, as the effective potential at zero temperature is a quartic polynomial irrespective of whether Ai — 0. It 
does however, affect the masses of the mesons. If Ai = 0, the following relation holds: ml — m^ = ^1^^ — m2 = 4M2, 
a property of the NJL model already noted in Ref. [s^. Clearly, the bosonized NJL model does not yield the most 
general linear sigma model. However, it gives additional contributions to the vacuum that usually are not taken into 
account in the linear sigma model coupled to quarks [22, 2§>] ■ In Ref. 29] it is noted that upon inclusion of fluctuations 
using an RG flow equation, the transition becomes second order. This boils down to including quark loop effects at 
zero temperature too and is consistent with our findings. 

To conclude, the mesonic part of this bosonized NJL Lagrangian is equal to the mesonic part of the LSMg model. 
So the mesons are treated in same way in the two models, but the vacuum contributions are treated differently. Since 
neither model is directly derived from QCD, it is not straightforward to draw a conclusion about the order of the 
phase transition expected in QCD. If the NJL model is viewed as a model for the microscopic theory underlying the 
low energy mesonic theory, it would not seem justified to neglect the logarithmic term at zero temperature. 

It is straightforward to bosonize the NJL model for 0^0 when also {^n^^ip) can become nonzero, leading to cross 
terms that mix the cr-field with the 77-field and OQ-field with the 7r-field, but we do not give the expressions here as 
they do not lead to any additional insights. 
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IV. CONCLUSIONS 



In this paper the high-T CP-restoring phase transition a,t 6 — n was discussed for two different models which aim to 
describe the low-energy QCD phenomenology, the NJL model and the linear sigma model coupled to quarks. Although 
the models are related, the philosophy of how the mesons are treated is quite different in both models. In the NJL 
model they are bosonic states of quark-antiquarks, whereas in the LSMg model they are the fundamental degrees 
of freedom, interacting with quarks at nonzero temperature. Using the bosonification procedure of Eguchi, one can 
show that a bosonized NJL model gives a linear sigma model, in which mesons are treated in the same way as in the 
LSMg model. However, the vacuum contributions arising from the quark degrees of freedom are different. The LSMg 
model was motivated for high temperatures, when constituent quarks are light and mesons are heavy. Therefore, it 
is assumed that quarks only play a role at nonzero temperature and do not affect the vacuum contributions at zero 
temperature. On the other hand, in the NJL model contributions by the quarks are necessarily taken into account 
also at zero temperature. The temperature dependent contributions to the effective potential are equal in both 
models, coming exclusively from the quarks. In the end, the effective potentials of the models only differ in their zero 
temperature contributions. Nevertheless, this directly affects the nature of the phase transition at high temperature 
a.t 9 — TT. 

The temperature dependence of the ground state of both models was investigated using a Landau-Ginzburg analysis. 
The difference between the models is that the potential as a function of the order parameter of the LSMg model 
contains a non-analytic logarithmic term, whereas the potential of the NJL model is a quartic polynomial near the 
phase transition. It is this logarithm that makes the difference, it affects the order of the phase transition. This 
logarithm comes from the contribution of the quarks at zero temperatures, but neglecting these contributions will 
affect the high temperature results qualitatively a,t 9 — t:. A similar effect occurs for the chiral symmetry restoration 
phase transition at 9 = close to the chiral limit, i.e. for sufficiently small explicit symmetry breaking. The absence 
of explicit CP violation is therefore an important aspect of the physics at = tt. 
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